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Double Asymptotic Scaling '96 * 
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We review recent HERA data on the structure function Fi at small x and large Q 2 . We show that the salient 
features of the data are revealed by comparing them to the double asymptotic scaling behaviour which F2 is 
predicted to satisfy in perturbative QCD. 



Double asymptotic scaling of the structure 
function F2(x,Q 2 ) in the two variables a = v^C 
and p ee TeTC (where f = ^ and C = ^g^) § 
is a simple consequence 0] of the QCD evolu- 
tion equations at small x and large Q 2 , and man- 
ifests the basic structure of asymptotic freedom in 
perturbative QCD ||. Because double scaling is 
spoiled if the dominant behaviour of F2 at large a 
and p is not governed by perturbative evolution, 
plots of F2 as a function of the scaling variables 
a and p were proposed in ref. 0] as a means of 
testing graphically the validity of the QCD pre- 
diction. Indeed, this analysis showed at a very 
early stage that, contrary to widespread expec- 
tation, perturbative QCD in the leading log Q 2 
approximation accurately describes the behaviour 
of the structure functions measured at HERA. 



While the availability of much more precise 
data on F2 now justifies a full NLO anal- 

ysis [|)-D based on the Altarelli-Parisi equations, 
whose applicability in this regime is now estab- 
lished beyond reasonable doubt, many instructive 
features of the small x, large Q 2 behaviour of F2 
are apparent from its double scaling behaviour. 

The basic double scaling behaviour of the data 
is displayed in Fig. 1 , which shows all data on F2 
which pass the cuts p, a > 1, Q 2 > 5 GeV 2 (so 
that the leading-order approximation is valid), as 
determined in the 1994 HERA run by HI j§ (a) 
and ZEUS (5) (b), plotted versus a (with xq = 0.1 
and Qo = 1 GeV), after rescaling of the leading 
one-loop sub-asymptotic behaviour [Q. Indepen- 
dence of In F2 from p (p scaling) is demonstrated 
by the fact that all data lie on the same line. 
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Fig. 2b 

In order to test for a scaling, however, we 
must make sure that the slope 2j of this line 
is that predicted in perturbative QCD, i.e. that 
the observed scaling is not an accidental conse- 
quence of the particular kinematic region covered 
by the HERA experiments. Note that the value 
of 7 = y/l2/f] = 6/^33 — 2nj is an absolute, 
parameter-free prediction of perturbative QCD. 
A simple way of doing so is to scale away the full 
asymptotic prediction, i.e. to rescale the data by 
the factor 

R F = Np^e-^+^'p, (1) 

with N an arbitrary normalisation. 

The result is displayed in fig. 2, where, as a 
cross-check, the dependence on p is also shown 
explicitly. Clearly both the HI (a) and ZEUS (b) 
data are in good agreement with the QCD predic- 
tion, which is attained asymptotically at large a. 
The (rescaled) data are now plotted on a linear 
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Fig. 3b 

scale: as a consequence, because of the definition 
of the scaling variables, deviations from the pre- 
dicted behaviour which are linear in x or log Q 2 
would appear as exponential scaling violations. 

For example, we can immediately exclude that 
the data behave as a scale-independent fixed 
power: in fig. 3 the HI data are compared to 
the function Ax~ x (A is an arbitrary constant) 
with (a) A = 0.08 and (b) A = 0.35.0 This is of 
some theoretical relevance because a (scale inde- 
pendent) power-like behaviour F% ~ would 
correspond to the dominance of a fixed Regge 
pole. The value A = 0.08 is that corresponding to 
the "soft pomeron" which dominates high-energy 
elastic proton-proton scattering g|, while A > 
0.35 is supposed to be [|l0| the smallest value one 

3 In this and subsequent plots the three curves correspond 
to p=l.l, 2.2, 3.3 on the a plot, and a = 1, 1.5, 2 on the 
p plot (dot-dashed, solid and dashed curves, respectively), 
which bracket the HERA kinematic region. 
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may obtain from a computation of the pomcron 
intercept based on the Lipatov equation in 
the presence of infrared cutoffs. Larger values of 
A would lead to yet stronger violations of double 
scaling. 

Of course such a power-like behaviour might in- 
stead only apply at some starting scale Qq as an 
initial condition to perturbative evolution. The 
leading asymptotic behaviour of F 2 then depends 
on the kinematic region Q: when p < ? the 
double scaling rise dominates over the boundary 
condition and F 2 displays double scaling, while 
for larger values of p the boundary condition is 
reproduced, up to calculable subleading correc- 
tions, i.e. [Q 

When A = 0.35 the transition occurs around p w 
3.5: it follows that fig. 3b allows us to exclude any 
boundary condition which is more singular than 
A « 0.35. Thus, the observation of double scaling 
excludes strongly rising boundary conditions |l]] . 

The behaviour eq. |2| corresponds essentially 
to approximating the Mellin-space anomalous di- 
mensions j lJ (N) which characterise the evolu- 
tion of singlet parton distributions with the value 
which they take at N = A. It is sometimes 
claimed that a fixed Regge pole might dom- 
inate the anomalous dimension (rather than the 
structure function F 2 , as assumed in producing 
the curves of fig. 3). The behaviour eq. || would 
then be applicable everywhere: of course, this 
would mean a breakdown of usual perturbative 



evolution in the region p < j where instead dou- 
ble scaling should hold. Comparison in fig. 4 of 
the resulting behaviour (eq. 2 with A = 0.35) to 
the data shows that this scenario is in any case 
clearly excluded by the data essentially because 
of its large deviation from the perturbative pre- 
diction as p decreases, and because the growth 
with a is too steep. 

Given that the data seem to exclude most of the 
commonly proposed alternatives to double scal- 
ing, one may wonder whether it is possible to 
mimic the perturbative double scaling prediction 
of eq. [l] by a simpler functional form. In ref. |l3| 
a good fit to the HI data B was obtained by 
assuming F 2 — a + b log f& log 21 , with a and b 

^0 x 

free parameters. The corresponding best-fit |l3[ 
curves are displayed in fig. 5, which shows that 
indeed there is reasonable agreement with double 
scaling in the HERA region. The agreement is 
only local, however: for example at large a and 
moderate p, i.e. at very large Q 2 and small x 
(where no data are available) the linear growth in 
In Q 2 leads to a discrepancy with double scaling, 
which corresponds instead to a rise of F 2 weaker 
than any power of InQ 2 . This parametrisation 
thus lacks any predictive power: the slope of the 
rise is fitted to the data rather than being pre- 
dicted, and thus cannot be reliably extrapolated 
outside the data region. It does however show 
that, due to the slow variation of £ and £ with 
Q 2 and x, the double scaling behaviour of eq. |l| 
admits simple local approximations. 

Despite its success in describing the data, it is 
clear from fig. 2 that double scaling is only correct 
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asymptotically. Indeed, a rise is visible in the p 
plot at large p, where Q 2 is low and the data get 
close to the boundary condition, while the data 
on the a plot appear to drop somewhat, except at 
the largest a values, indicating that the slope of 
the rise of F2 is sub-asymptotically smaller than 
the predicted one. 

We may investigate whether these violations 
of double scaling may be understood in a sim- 
ple way by seeing whether they are reduced by 
the inclusion of a two loop correction to the lead- 
ing double scaling behaviour. The correction may 
be computed (w) by determining the NLO cor- 
rection to the leading singularity of the anoma- 
lous dimensions in N (which has the same loca- 
tion and strength as at leading order), and de- 
riving the corresponding two loop asymptotic be- 
haviour, which turns out to be scheme indepen- 
dent. We can then once again test graphically 
this two loop correction by rescaling the data to 
the two loop asymptotic prediction, i.e. with the 
rescaling factor 
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where (3i is the two loop coefficient of the /3 func- 
tion, and Rf is given by eq. [I], with £ computed 
using the two-loop expression for a s . Fig. 6 shows 
that indeed the two loop correction eq. |J removes 
the main systematic violation of double scaling, 
namely the sub- asymptotic fall in the a plot. The 
two loop correction also has the effect of mov- 



ing upwards the optimal reference scale Qo, to 
Qo = 1.8 GeV (this is the main cause of the 
change of scale in the plots of fig. 6 compared 
to the previous one-loop plots). 

However, the overall p scaling deteriorates 
somewhat. Even though part of the p scaling vi- 
olations, in particular those at large p, are pre- 
sumably the result of boundary effects, this sug- 
gests that the part of the two loop anomalous di- 
mension not included in the determination of the 
asymptotic behaviour eq. |] may actually be im- 
portant in most of the HERA region. Otherwise, 
stated, this suggests that the summation of logs of 
Q 2 may still have a leading role over the summa- 
tion of logs of i . This is consistent with the result 
of ref . (?]] , where it was shown that the observed 
perturbative evolution of F2 is most accurately 
described by two-loop Altarelli-Parisi evolution, 
while the inclusion of the summation of logs of 
- in the evolution equation inevitably makes the 
agreement with the data worse. 

This also suggests that a more accurate descrip- 
tion of the data requires a full two loop computa- 
tion. The results of such a computation [0J|] are 
also displayed in fig. 6. The curves here corre- 
spond to p—2, 3, 4 on the a plot, and <r = 1, 1.2, 
1.4 on the p plot (dot-dashed, solid and dashed 
curves, respectively). The remarkable accuracy of 
the two loop computation is however best shown 
by plotting the ratio R of the data to the best fit 
versus the two scaling variables (fig. 7) . This fit is 
obtained with only two free parameters, namely 
the exponents X q and X g which characterise the 
small x behaviour x~ Xi of the input singlet quark 
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Fig. 8 

and gluon distributions (see ref. || for details of 
the fitting procedure). This is possible because 
the dominant double scaling behaviour is largely 
independent of the details of the input parton dis- 
tributions jjj. In fact, in most of the HERA re- 
gion any starting distribution leads to essentially 
the same behaviour, provided only that A^ are not 
too large (i.e. A < 0.35, as discussed above). 

The availability of data at low x and Q 2 (i.e 
large p and large a) however, allows now to put 
more stringent bounds on the acceptable values 
of the exponents \ than the asymptotic double 
scaling analysis alone would allow. Indeed, the 
full two loop fit H leads to the values \ — 0.25± 
0.02 and \ g = 0.09 ± 0.09 at Q = 2GeV in the 
MS scheme: hence the values, say, X q = X g = 
0.35 are excluded by several standard deviations, 
even though they would only lead to violations of 
double scaling at large p > 3.5. This is apparent 
in the plot of the data/theory ratio determined 
with these values (fig. 8): the excessive rise of 
the boundary condition leads the computed F2 
to overshoot the data close to the boundary, and 
thus to a systematic fall as p increases. 

The precise values taken by the exponents A.; of 
course depend on the particular parameterisation 
adopted, on the factorization scheme and on the 
value of the starting scale, and thus do not have 
any direct physical meaning. It is nevertheless 
significant that the data require that at low Q 2 
singlet parton distributions, and the gluon distri- 
bution in particular, do not rise strongly with -, 
in qualitative agreement with expectations based 
on the soft pomeron of Regge theory B . 



The description of scaling violations afforded 
by standard two-loop Altarelli-Parisi evolution is 
now accurate enough that scaling violations due 
to the summations of logs of i may be excluded 
phenomenologically |j. This, together with the 
success of double scaling, which implies the ab- 
sence of strong nonperturbative boundary effects 
which would spoil the perturbatively generated 
behaviour, opens the way to very accurate QCD 
phenomenology at HERA, such as, for instance, 
precision determinations of a s A complete 
theoretical understanding of the underlying rea- 
sons for this happy state of affairs is however still 
missing, and will certainly require going beyond 
the simple picture of double asymptotic scaling 
summarised here. 
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